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Spin motive force is a spin-dependent force on conduction electrons induced by magnetization 
dynamics. In order to examine its effects on magnetization dynamics, it is indispensable to take 
into account spin accumulation, spin diffusion, and spin-flip scattering since the spin motive force is 
in general nonuniform. We examine the effects of all these on the way the spin motive force generates 
the charge and spin currents in conventional situations, where the conduction electron spin relaxation 
dynamics is much faster than the magnetization dynamics. When the spin-dependent electric field 
is spatially localized, which is common in experimental situations, we find that the conservative part 
of the spin motive force is unable to generate the charge current due to the cancelation effect by the 
diffusion current. We also find that the spin current is a nonlocal function of the spin motive force 
and can be effectively expressed in terms of nonlocal Gilbert damping tensor. It turns out that any 
spin independent potential such as Coulomb potential does not affect our principal results. At the 
last part of this paper, we apply our theory to current-induced domain wall motion. 
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In a ferromagnetic system, dynamics of space-time de- 
pendent magnetization vector M(x, t) is described by the 
following phenomenological equation [H-Q 
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This is the Landau-Lifshitz-Gilbert (LLG) equation gen- 
eralized to include the spin-transfer torque terms (last 
two terms) . Here, H e f f is the functional derivative of en- 
ergy with respect to M, a is Gilbert damping constant, 
M s is saturation magnetization, (Ib is Bohr magneton, 
e is electron charge, is nonadiabaticity, and j s is spin 
polarized electric current given by difference between cur- 
rent by spin-up and spin-down electrons. LLG equation 
describes the dynamics of magnetization under applied 
electromagnetic fields. On the other hand, there exists 
a reciprocal process; temporal and spatial variation of 
magnetization induces additional electromagnetic fields 
on conduction electrons. These fields are spin-dependent 
and in general nonconservative. Thus the resulting spin- 
dependent motive force is called spin motive force (SMF). 
SMF is firstly predicted by Berger [5[ and recently for- 
mulated by generalizing Faraday's law [6]. It is also sug- 
gested _p that SMF is also described by spin pumping 
effect M. Recently, SMF and its effect is intensively stud- 
ied 0,H-i3 in this field. 

Without any other perturbation, the explicit expres- 
sions of the induced spin electromagnetic fields are known 



where f an d 4 stand for spin-up and down electrons. 
The fields generate spin-dependent Lorentz force FJ^ — 
e^J^ + vU x B|^). While the term "motive force" refers 
to quantities of voltage dimension, the term SMF is some- 
times used to denote FJ^. In this paper, we adopt the 
latter terminology. It is easily noticed that SMF is spin 
dependent, nonconservative, spatially varying, and local- 
ized in general situations. As one can see from Eqs. ([2]) 
and ©, SMF is usually too small to be measured di- 
rectly. Another way is to study the effect of SMF on 
magnetization dynamics. Since SMF induces additional 
spin current, additional spin-transfer torque arises and 
it changes LLG equation. Consequently, LLG equation 
(without applied spin current) is modified as Q 
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where T> is modified damping tensor given by 



T>u = ad,. 
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and r\ = [iB^a /2e 2 M s . Here, a is electrical conductivity. 

However, the previous work has a crucial limitation 
that spin density has been considered to be constant 
while spin density in reality is nonuniform since SMF is 
spatially varying. The main consequence of the nonuni- 
form spin density n^(x, t) is diffusion current propor- 
tional to -Vn^, which suppresses the effect of SMF. 
Therefore, a realistic model should take into account 
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spin accumulation (nonuniform spin density), diffusion 
and spin-flip scattering. The purpose of this paper is to 
find the spin density rv^, diffusion current, total current 
induced by SMF, and their effect on magnetization dy- 
namics in the presence of spin accumulation, diffusion, 
and spin-flip scattering. As one shall see in Sec. IIII1 the 
solution of rv^ in the most general situation is too com- 
plicated to study the effect on magnetization dynamics. 
To obtain simple analytic expressions, we take an ap- 
proximation that spin-flip time is much shorter than the 
time scale of magnetization dynamics. As a final com- 
ment, our result does not assume any specific form of 
SMF. Thus, it remains valid for the modified SMF due 
to, for instance, nonadiabaticity (T^j, spin-orbit coupling 
[201 ] , and other kinds of spin dependent electric field [22| . 

Several previous works are closely related to our work. 
Spin drift-diffusion equation, which has similar form to 
our theory is suggested in Ref. And, the effect on 

spin and charge current is investigated from Boltzmann 
equation in Ref. We set our starting point as the 

equation of motion of conduction electrons in Ref. [l| 
to make our analysis consistent with previous theories in 
this field. Different from the previous theories focusing 
on ID, we successfully generalized our result to 3D, and 
found that nonconservative part of SMF plays a crucial 
role in current in a higher dimensional system. In Sec. 
IIII1 we compare our result with the previous theory qual- 
itatively and quantitatively In addition, we investigated 
the effect of charge neutrality on our results. It turns 
out that charge neutrality potential does not change our 
principal results, charge current and spin current, even 
though it changes charge density and spin density. Fur- 
thermore, we show that any spin independent potential 
cannot alter our principal results, either. 

This paper is organized as follows. In Sec. UH we 
construct the spin drift-diffusion equation and introduce 
variables. Then, we solve the equation in Sec. IIII| and 
discuss various implications. In Sec. IIV| we apply our 
result to current-induced domain wall (DW) motion and 
briefly discuss the effect of spin diffusion. In Sec. [V] we 
generalize our theory for general boundary condition and 
for general spin indendendent potentials. Finally, in Sec. 
IVI1 there are concluding remarks. Technical details are 
in Appendices. 



II. MODEL 

A. Spin drift-diffusion equation 

To construct the equation of n^(x,t), we take the 
starting point as the equation of spin density m of con- 
duction electrons [l|, 
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Here, J is spin current tensor, r ex = h/SJ ex , and S de- 
notes the magnitude of spin of local magnetization. The 



left-hand side is based on the continuity equation. The 
first term on the right-hand side is the precession term 
due to the exchange coupling between conduction elec- 
trons and magnetization. (T) includes the effect of spin 
scattering processes. Here, the second rank tensor J is 
defined by 
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The effect of the perpendicular component to M of Eq. 
([6|) is already investigated by Zhang and Li |1[ , and they 
found the nonadiabatic term of LLG equation. In the 
absence of spin accumulation, the magnitude of m is 
constant, so it suffices to solve only the perpendicular 
component of the equation. However, in the presence 
of the spin accumulation, the magnitude variation of m 
should be also studied. We define spin number density 
n s = \m\/ fiB- Taking care of the fact that n s has space- 
time dependence, the parallel component of Eq. ([5]) to 
rh = m/|m| results in 
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It is convenient to separate the variables to that of up 
and down electrons. n s = n' — rv' and j s = y —y. Here, 
and denote spin number density of spin- up/down 
electrons and charge current density generated by spin- 
up/down electrons, respectively. Equation ((SJ) is nothing 
but the continuity equation containing spin nonconserv- 
ing processes described by T. To obtain independent 
equations of spin- up/down electrons, we use the follow- 
ing continuity equation of total electron number density 

dn e 1_ . . . 

-^ = -V-j e , 9 
Ot e 

where n e — rv + is electron number density and j e = 
jt _|_ ji j s charge current density. Combining Eqs. and 
©, one obtains 
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Note that T represents spin-flip scattering processes. As 
a simple model, we take the well-known form of spin-flip 
scattering, 



1 TV rr 
rh ■ (T) = —r -. 
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where is characteristic time of the spin-flip scattering 
process from spin-up to -down state, and is similarly 
defined. Then, 

which is the spin drift-diffusion equation. Similar form 
of Eq. (fT2"]l was also suggested in Ref. [l(| • 
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For simplicity, we may assume without losing general- 
ity that the SMF is turned on at t = and that, for t < 0, 
the system is in equilibrium. We set n^(x, t = 0) = rv Q , 
where nj^ is equilibrium electron density of spin up and 
down at t < 0. By definition, the equilibrium density 
nj^ is the equilibrium solution of Eq. (12) for t < 0. 
Inserting n^x, t < 0) = to Eq. (12), one ob- 
tains an important constraint tIq/t^ = Uq/t^. With 

i* i 

the help of this constraint, four variables n , Hq, t-j- and 
"H - can be described by three variables, nj, tt-q, and t s j 
(t^ 1 = t*~ 1 + r i_1 ). Then, Eq. (Til]) is rewritten with 
only one spin- flip time r s f . In addition, current can 
be written as a'^'E'l^+eD'^'Vn'^, where and are 
respectively the conductivity and SMF (divided by e) for 
spin- up/down electrons. Then, one straightforwardly ob- 
tains the final form of the equation of our model. 



Variable definitions and relations 



In Sec. IIIIl there appear several variables and quan- 
tities which have not been defined yet. To help readers, 
we present definitions of them here, rather than Sec. IIIIl 

Since Eq. (fl3l) is coupled, it is convenient to solve it in 
matrix form. Hence, we define spin accumulation vector, 
which is a column vector defined by 



Similarly, we define current density and SMF vector. 
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As mentioned in Sec. [fl we treat E^ as nonconserva- 
tive, spatially varying fields. In addition, it is assumed 
that spin dependence of E s is given by E| = — E^ = E s . 
Slight generalization of our theory at the final step al- 
lows to investigate the formula for E| ^ — E£. No other 
restriction of Ej* is not assumed in order to obtain max- 
imally generalized result. 

As suggested in Ref. 10], in realistic systems, 
the Coulomb interaction should be taken into account. 
Hence, one introduces Coulomb potential V c and add it to 
the spin motive force as E|^ — > E|^ — W c . The Coulomb 
interaction strongly suppresses the charge accumulation. 
Mathematically the interaction may thus be handled by 
imposing the charge neutrality constraint. We show in 
Sec. |V]that charge neutrality constraint changes electron 
densities, but not currents. Hence, the LLC equation is 
hardly affected by charge neutrality potential. For this 
reason, we do not take into account the Coulomb inter- 
action until Sec. |V]in order to show simple logical flow. 

Note that all variables in Eq. (|13l) are not indepen- 
dent. Einstein's relation is given by a^l = e 2 Z)tl]ytl 
where ]SP+ is density of states of spin- up/down elec- 
trons at Fermi energy. Since cx n^, one obtains 
,jt /D^tIq = / D^iiq. This is one of the key constraints 
of our model. 



The solution of Eq. (|T3|) is very complicated as one 
shall see in Sec. IIIIl To gain an insight, it is illustrative 
to assume that r s f is much smaller than the time scale 
of magnetization dynamics so E s is almost constant in 
time scale within T s f. We found that, in this limit, the 
solution is much simpler and it is easier to catch physical 
meanings. 



Equations pT| - (jT6|) are related by the following relation. 

The first term in right-hand side corresponds to conven- 
tional electrical current and the second term corresponds 
to diffusion current. 

Instead of diffusion constants, it is more physical and 
intuitive to express results in terms of spin-flip length 
which is defined by X^j = yj D^^r s f . The averaged spin 
diffusion length is also defined by the conventional way 
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By Einstein's relation, Eq. (fT8|) is equivalent to 
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Combining Eqs. (TT5]) and (fT5|) . A^i is represented in 
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where a = + is total electrical conductivity. 

Conductivity polarization P and density polarization 
P n are defined by 



P = 

Pn = 
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With these polarizations, er^ and nj^ are represented 
in terms of a and (nj + Uq) as er^ = (1 ± P)cr/2 and 
n^ = (l±P n )(nl + 4)/2. 
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Lastly, we use a mathematical convention that A(k) is 
the Fourier transform of a position dependent function 
.A(x) with respect to x. That is, 



i(k) ^L4(x)](k) 



1 
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for a (i-dimensional system. 
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III. CHARGE AND SPIN CURRENTS IN THE 
PRESENCE OF SPIN DIFFUSION 



A. Solution of the spin drift-diffusion equation for 
localized electric field 



Before solving Eq. (JT3J) for general cases, we first solve 
the equation for localized E s since E s is localized in most 
cases. In Sec. |Vj we generalize our theory to include 
spatially extended E s . 

Since E s is localized, it is possible to take Fourier trans- 
form with respect to position. Then, £ and Af are well- 
defined localized functions of k except initial condition 
part. In addition, localized E s implies that the bound- 

because 



,T4- 



ary condition is given by n^(|x| — > oo,i) 
E s (|x| — > oo,i) = does not affect spin density. After 
Fourier transform, Eq. (|13p is written as, in matrix form, 



~dt 
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n = 



(24) 
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Equation (|24l) is a first order ordinary differential equa- 
tion with respect to t and the initial condition is given 
by Af(k,t = 0) = (2Tr) d / 2 S(k)(nl n£) T . The solution is 
simply given by 



-nt/r Bf 




■AfCM) 



Since the first term of Eq. (|26|) represents the time vari- 
ation of equilibrium number density, one can realize that 
the term should be given by e - Qt / T 'fJV(k, 0) = A"(k, 0). 
Mathematical derivation of this argument is given in Ap- 
pendix [3] The second term of Eq. (121)1) is almost impos- 
sible to take inverse Fourier transform. Hence, as men- 
tioned, we use an approximation that r s f is very small. 
In this limit, Appendix IB1 shows that 
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where Q(t) is Heaviside step function. By this approxi- 
mation, solution of the spin drift-diffusion equation Eq. 



(1261) becomes 
Sf(k,t) =A/"(k,0) 
The inverse of is explicitly given by 
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Now, excited charge density An e = (rfi + n^) — (nj + Uq) 
and excited spin density An s = (rft — n}) — (nj — Uq) are 
given by, 
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in k-space. 

At this stage, there is no need to show complicated real 
space expressions of the densities, because what affects 
LLG equation mainly is spin current. In the next subsec- 
tion, we find the expressions of charge and spin currents 
in bith fc-space and real space. 



B. Charge and spin currents 

Charge and spin currents in the absence of spin diffu- 
sion are given by 



c^El + cr^Et = PaE s 
c^Et - cr+Et = ctE,. 



(33) 
(34) 



In this subsection, how the spin current and charge cur- 
rent generated by SMF is changed by spin diffusion from 
Eqs. (|3"3")l and (|34l) is examined with the help of Eqs. 
(TiTl) and ([281). After some algebra, 
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Now, the expressions of charge current j e 



spin current j 



■ y are straightforward. 
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We present j s in terms of j e in order for one to see easily 
js = je for P = ±1 ; for perfectly polarized electrons 
without spin-flip, the spin current should be the same 
amount of the charge current. 

Equations (j3"oT) and (|3"7]l are the principal results of this 
paper. One can obtain e?-dimensional real space expres- 
sions by taking inverse Fourier transform. This is one 
of the key advantages of our theory. Our result is easily 
generalizable to rf-dimcnsional result. As one shall see, 
it turns out that the nonconservative part of E s plays a 
crucial role in a higher dimensional system. 

First of all, we present ID result. ID real space ex- 
pression of Eqs. (|3"6")l and (13"T1) are 



j e {x,t) = 0, 



(38) 



j s (x,t) = (1-P 2 ) / dx'— — aE s (x',t). {39) 

One can notice that the ID charge current is perfectly 
canceled by diffusion current for small spin-flip time 
limit. This is natural in the sense that, for small spin- 
flip time, the system tends to be in equilibrium at each 
time t. At equilibrium, the gradient of chemical potential 
vanishes, so does charge current. However, spin current 
does not vanish by this reason because of spin noncon- 
serving process. Due to spin diffusion, it is natural that 
the spin current becomes nonlocal with integration kernel 
width A s /. Here, the factor (1 — P 2 ) seems unexpected. 
This factor comes from spin diffusion effect, and should 
exist regardless of diffusion length scale. It yields more 
cancelation for more polarized electrons. Eventually, for 
P = ±1, spin current also vanishes, which is actually 
required since j s = j e in the limit P = ±1. 

Equation (|39p behaves quite differently for two limiting 
cases. Let A be the characteristic length scale (such as 
DW width) of localized E„. If A s/ < A, e l*-*'l/-W « 
2X sf S(x-x'). Then, 



j s (x,t)^(l-P 2 )aE s (x,t), 



(40) 



which is local. Very short diffusion length cannot make 
the spin current nonlocal. It is very interesting that (1 — 
P 2 ) factor does not disappear even though diffusion effect 
is very small. The existence of diffusion makes (1 — P 2 ) 
factor regardless of how the effect is strong or weak. For 
Kf > A, 
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FIG. 1: (color online) Results of micromagnetic simulation 
for (a) charge current and (b) spin current induced by SMF 
generated from ID DW oscillator, which is fixed at x = 0, but 
rotating with to = 10GHz, in the presence of spin diffusion 
with A s / = 0.5nm (black rectangle), 5.0nm (red circle), and 
50nm (blue triangle). Solid line represents currents without 
diffusion. Here, DW width is set to be lOnm. 



where E s is averaged SMF and X(t) is the position of 
localized E s (such as DW position). One can see that 
the current is highly suppressed by the factor A/A s /. In 
this highly diffusive regime, spin current is also highly 
suppressed. 

The main features of our result is similar to those of 
Ref. except for vanishing charge current. They 

claim that charge current can exist in general, while we 
find that Einstein's relation prevents the existence of 
charge current in ID. 

We confirmed our analytic result by comparing it with 
micromagnetic simulation for a ID DW oscillator. The 
result is in Fig. [TJ General features are the same as an- 
alytic result ; i) charge and spin currents are highly sup- 
pressed by spin diffusion, ii) spin current becomes non- 
local, iii) charge current almost vanish independently of 
diffusion length, and iv) spin current is more suppressed 
by larger diffusion length. 

Now, we generalize the results to higher dimension. In 
2D and 3D real spaces, Eqs. (|36| and (j37l are converted 
to 
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respectively. Here, -Ko^) is the zeroth order modified 
Bessel function of the second kind. One can be uncom- 
fortable because j 2D seems dependent on A s / while Eq. 
(|36[) is not. Here A s / is included because the argument 
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of logarithmic function should be dimensionless. In fact, 
one can easily check that jI d does not depend on A s / by 
using In |x — x!\/\' s f — In |x — x'|/A s /+(constant) for any 
positive X' s f. 

From now on, we present only 3D expressions but omit 
2D, for simplicity. One can easily obtain 2D expressions 
by replacing the integral kernels l/|x — x'| — > — 2 In |x — 

x'|/A s/ and exp(-|x-x'|/A s/ )/|x-x'| -> 2K Q {^j 1 )- 

The reason why the expressions depend on dimension 
is nothing but the fact that the inverse Fourier trans- 
forms, which give integration kernels, depend on dimen- 
sion. Hence, essential physics are the same for 2D and 
3D except the mathematical expressions of the integra- 
tion kernels. 

Overall features of the spin current in higher dimension 
is similar to ID, except for the existence of nonvanishing 
Pj e term. However, the main feature of j e is completely 
different from ID case. First of all, j e does not vanish. 
We argued qualitatively why the charge current vanishes 
in a ID system by using chemical potential argument. 
However, in a higher dimensional system, chemical po- 
tential cannot be defined in general since E s is noncon- 
servative. Note that diffusion current D^Vn'^ is con- 
servative. Note also that nonconservative field cannot be 
canceled by conservative field. This is why the charge 
current exists in a higher dimensional system. The non- 
local term in Eq. (l4"2l can be interpreted as Coulomb 
potential under charge density —V • E s . The canceled 
part of the charge current is nothing but conservative 
Coulomb part of E s . Secondly, it is very interesting that 
Eq. (|4"2"j) is converted after some algebra to 



(a) 



Pa 

47T 



^ V x (V x E, 



(44) 



Now one can notice the importance of nonconservative 
part of E s (nonvanishing V x E s ) and the dependence 
of this on the charge current. If V x E s happens to be 
zero, the charge current also vanishes, and this is consis- 
tent to the chemical potential argument. Lastly, it is also 
interesting that charge current does not depend on dif- 
fusion length. This is qualitatively understandable from 
the fact that the effect of E s is maximally canceled by 
diffusion current in small spin-flip time regime, regard- 
less of spin diffusion length. Here, maximal cancelation 
is slightly different from perfect cancelation in ID case in 
the sense that the (conservative) diffusion current cannot 
cancel PcrE s perfectly in principle. 

It might be ambiguous what the "conservative part" 
of a vector field is mathematically. Helmholtz's theorem 
guarantees that a spatially localized vector field can be 
uniquely decomposed into conservative (curl-free) part 
and solenoidal (divergence- free) part. Note that the sec- 
ond term in Eq. (|4"2")l is exactly the same as the formula of 
(negative of) conservative part of the Helmholtz decom- 
position. Note also that the resulting total current [Eq. 
(|44p] is divergence-free. Therefore, the charge diffusion 
current and total charge current are respectively given 
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FIG. 2: (color online) Results of the micromagnetic simu- 
lation for vortex resonant oscillation in 2D thin film. Here, 
the vortex core size is about 5nm, the resonant frequency is 
605.5MHz, and P — P n = 0.7. The size of arrows are log- 
scaled (size oc ln(l+norm)). (a) Magnetization profile at the 
time when all calculation was performed. The blue arrow 
denotes the direction of vortex core motion, (b) SMF. (c) 
Accumulated charge density — eAn e (divided by <rr s //A^). 
The maximum value (white) is 1.39 x 10 _5 V and the mini- 
mum value (black) is -1.62 x I0 _5 V. (d) V ■ E s . The max- 
imum value (white) is 8.43 x 10 10 V/m 2 and the minimum 
value (black) is —7.42 x I0 10 V/m 2 . (e) Charge diffusion cur- 
rent, which is the conservative part of SMF. (f) Total charge 
current, which is the solenoidal part of SMF. 



by conservative part and solenoidal part of Helmholtz 
decomposition. One shall see in Sec. |V]that this claim is 
generally valid for arbitrary boundary conditions. 

To understand the maximal cancelation of the charge 
current qualitatively, it would be very helpful to visualize 
the Helmholtz decomposition of the charge current. We 
performed a micromagnetic simulation for vortex reso- 
nant oscillation in 2D thin film. When the vortex DW 
wall [Fig. Ufa)] moves along the blue arrow, SMF is gen- 
erated as shown in Fig. [2j[b) . The spatial dependence of 
SMF induces charge accumulation as shown in Fig. Etc) . 
Due to the charge accumulation, charge diffusion current 
is generated as shown in Fig. ^[e) . Figure ^e) can be 
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qualitatively understood by Fig. [2fd). Recall that the 
diffusion current is given by Coulomb field generated by 
charge density —V • E s . Thus, dipole-like nature of the 
charge density — V-E s [Fig. [2fd)] implies dipole-like field 
[Fig. EJe)]. Summing up Figs. [2fb) and[2je), one obtains 
total current as Fig. [2If)- Note that the total current is 
definitely nonconservative (has finite curl) . It is interest- 
ing that Fig. [5Jf ) is similar to magnetic field generated by 
two separate conducting wires. This infers the solenoidal 
nature of the total charge current. It is very interesting 
that total charge current behaves as magnetic field rather 
than electric field. 



C. Magnetization dynamics and nonlocal damping 
tensor 

It is also important to see how the magnetization dy- 
namics is changed by our result. By analogue of Ref. 
Q, it is obvious that the modified LLG equation is de- 
scribed by Eq. ([T]), in which Eq. (|43| are added to j s 
terms, rather than Eq. However, there exist two 

nontrivial features in the modified LLG equation. 

The first one is spatial dependence of 7 and a. It is 
important to notice that 7 and a in Eq. ([1} are renor- 
malized parameters [l[, 



7o 1 + 



7o , 
— OiO 

1 



n s ^B 



1 



M s 1 + /3 2 
,n s /J,B 1 



M s 1 + /3 5 



(45) 



(46) 



where 70 and ao are original parameters of the system. 
There has not been any problem of this renormalization 
without spin accumulation, but n s is no longer constant 
in the presence of spin accumulation. Hence, 7 and a 
cannot be a simple constant in principle. At this stage, it 
is more convenient to write down LLG equation without 
parameter renormalization, 



1 



1 ti s ij,b\ 



1 + /3 2 M ( 

P 1 n s /j,B 
a l + f3 2 M s 



eMs 



(is ■ V)M 



P^b 
eM? 



= -70M x H e// 

Ms dt 
Mx (j S ' V)M. (47) 



In fact, the effect of n s variation is small. Note that 
n s fi B /M s = (nl-n^)fi B /M a x An s /(nJ-r4). Here, the 
first factor is of the order of 10~ 2 and the second factor 
is first order in SMF. Then, the effect should be very 
small compared to ordinary first order effect of SMF. In 
addition, one shall see in Sec. IIV1 that symmetry can 
reduce the effect of An s in collective coordinate level. In 
the case, the effect of An s vanishes through odd function 
integration. 

The next one is effective damping constant. Without 
applied electric field, the modified LLG equation is ob- 
tained by taking Gilbert damping as a damping tensor 



like Eq. (HJ) . Then, it is interesting to see how the damp- 
ing tensor is generalized by spin diffusion effect. In the 
presence of spin diffusion, local damping tensor T>ij be- 
comes nonlocal damping tensor 2?jj(x, x'). By using Eq. 
(|43|). it turns out that the damping tensor in Eq. (UJ is 
modified as 



T>ij(x, x') = aSijS(x — x ) + 



,, (1 -P 2 ) r?e -N- a: 'IA./ 



Ml 2\ sf 
x(Mx3 I M) i (M'x9 I M') j , (48) 



for ID and 

X>jj(x,x) = a5ijS(it — x 



n 



■£jm(x,x') 



aMj 

x(M x d k M)i(M' x diM.%, (49) 

P 2 a 1 
5y(x,x) = cr<Jytf(x-x) + -j — didj-. - 

47T X — X 



47T 



didj- 



X - X' 



(50) 



for 3D. Here, M' = M(x',t) and • is redefined as inner 
product with respect to both coordinate basis and spatial 
basis, 

[V ■ f(x)] 4 = J2j dVPtftx.xO/^x'), (51) 



for a vector field f . As a passing remark, the tensor D is 
indeed a damping tensor in the sense that it decreases to- 
tal magnetic energy. It can be demonstrated by showing 
energy dissipation is negative, 



dE_ 
~dt 



I 



H e// — ),/', 0. ,.-,2! 



The mathematical details related to Eqs. (|48|) - (l521) arc 
in Appendix [U] 



IV. EXAMPLE : DOMAIN WALL MOTION 

In this section, we apply our theory to ID DW mo- 
tion. We find equation of motion of collective coordinates 
(X(t),(j)(t)) of tail-to-tail transverse DW. Here, X(t) is 
DW position and <fi(t) is tilting angle. Mathematical de- 
tails of obtaining the collective coordinate equation is 
described in Appendix [Dl 

Without SMF, DW motion is described by 0, 



d(j) adX 
~~dl + X~dT 
1 dX d(f> 



A dt 



dt 



A ' 
b°j iK d 
' A Ms 



sin 2(f), 



(53) 
(54) 



where b® = Pj^Lg/eMs for applied current j, K4 repre- 
sents dipole field integration, and A is DW width. In the 
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presence of SMF, Eqs. (|53j) and (|54|) are modified as [2 



dcj) a dX 
~~dt + X~dt 
IdX dcj> 
X~dt 



A 



2firi dcj) 



a 



dt 



3A 2 dt 
>° 2rj dcj) 
A~ ~ 3A 2 ~dt ~ ~M. 



(55) 
sin 20.(56) 



One can expect that the effect of 77 will be suppressed by 
spin diffusion, so a renormalized parameter r) will replace 
r\. Thus, the expected equations of motion are 



dt 
1 dX 
X~dt 



adX 
X~dt 

d<t> 



X 

b l 

X 



2 fir} dcj) 
~ 3A 2 ~dt ' 
2rj dcj) _ jK d 
3A 2 dt M s 



(57) 
sin 20.(58) 



As described in Appendix IDl it turns out that Eqs. (|57|) 
and (|58j) are indeed valid, and the renormalized SMF 
parameter is given by 



fj = (l-P 2 )F(Ov, 



r(z) 



(59) 
(60) 



where £ = X/X s f and T(z) is the Gamma function. 

The renormalization is largely dependent on relative 
magnitude of DW width and spin diffusion length. The 
asymptotic behavior of the function F is given by 



F(0 = 



c 



3/-2 
5? 



8 T 

+ 0(C 



0(C 4 ) 

4 ) 



C«i, 
C»i. 



(61) 



For A <C X s f, Eqs. (|55| and (|56|) are reproduced except 
(1 — P 2 ) factor which should exist as already discussed. 
For A » A s /, the effect is highly suppressed by spin dif- 
fusion, so the overall effect of SMF goes as tt— rather 



than tj, 



AA. 



FURTHER GENERALIZATIONS 



A. Extended electric field - general boundary 
condition 



In the case that the magnetization dynamics is gen- 
erated by applied electric field, the electric field is no 
longer localized. Hence, it is necessary to generalize our 
theory to non-localized electric field ; E s does not vanish 
at |x| — > 00. In this case, £ does not give well-defined 
Fourier transform, but includes delta function parts. Fur- 
thermore, in order to obtain Eq. (j44jl . one obtains ad- 
ditional boundary terms when integrating Eq. (|42[) by 
parts. Hence, the charge current may not be canceled by 
the diffusion current even if the electric field is conserva- 
tive. 

The problem can be treated by Green's function 
method with given boundary condition, as described in 



the last part of this section. However, it is hard to catch 
the physical meaning, so we present more intuitive analy- 
sis for this case. It is very convenient to use the linearity 
of our theory. Suppose can be decomposed into two 



jmponents E^ = Ej 



■E 



Then, since Eqs. (Tl3|) and 



(fT"7f are linear, the current can also be decomposed 
into two components = j| + j| , where jj^ is the 
current generated by Ej. 

To take advantage of this linearity, we decompose E 
into three components as follows; any vector field E 
can be uniquely decomposed into irrotational part Ei, 
solenoidal part E2, and boundary part E3, 



E n = E n 4 

V • E^ ^ 0, 



E 



n 



E 



Ti 



V-E'. 
V-E 



2 4 ' 

Ti 



x E^ = 0, 



(62) 



n-El J -| xe8 = 0, (63) 



V x Ef - ^ 0, n • E^| xea = 



V x E 



n 



0, n-E^| xea ^0. 



(64) 
(65) 



Here, d means boundary and n denotes a normal unit 
vector perpendicular to the boundary. From given E^, 
one can obtain E^ by solving Laplace equation with 
Neumann boundary condition and E|^ 2 by Helmholtz's 
theorem. As previously mentioned, our result [Eq. (|44p ] 
indicates that the conservative source part E^ cannot 
contribute to the charge current while the nonconserva- 
tive solenoidal part E|^ can give rise to a nonvanishing 
contribution o^El + cr+E^. [Eq. @2J) or Eq. {77} for 
general boundary condition] 

In the presence of nonvanishing boundary condition, it 
is important to investigate the effect of Ej^ to the spin 
accumulation. Fortunately, the source term of Eq. (fT3|) 
depends on the divergence of E|^ only. Hence, Eg^ can- 
not contribute to the spin accumulation since V-eJ 4 " = 0. 
Therefore, E| can only contribute to the currents via the 
first term in Eq. (|17p . Consequently, the expressions of 
charge and spin current should be cr^Eji cr^E^ added to 
Eqs. (O and (03J (or Eqs. ® and (32)). 

It can be a good example that a constant spin inde- 



pendent electric field E 
wire. In this case. 



E, 



-pTi 
^3 



E n 



spin current should be E app ± a^E app 
([581) and (|59"1). so 

j e = crE n 



a ppX is applied in a ID 
x. Hence, the charge and 
added to Eqs. 



3s 



J app> 



PcjE n 



(1-P 2 



dx'- 



-\x—x' |/A S / 



2A 



(66) 

-aE s (x'im 



In this case, the charge current is not canceled by diffu- 
sion current, but only charge current generated by SMF 
(which is spatially localized) is canceled. 



B. Charge neutrality and other spin independent 
potentials 

There is another important physical consequence of 
spin accumulation we have ignored. This is Coulomb 
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potential of accumulated electrons. Since Coulomb po- 
tential is in general strong, electron tends to make local 
charge neutral. Hence, giving a charge neutrality con- 
straint, 



An? + An 1 = 



(68) 



is a good approximation. To do this, it is convenient 
to introduce charge neutrality potential Vr and replace 
E ti _> E |i _ vv c as suggested in Ref. Then, the 

solution of the spin drift-diffusion equation [Eq. ([25)1 ] is 
modified by 



AfiT 
An 1 



e { a± J * k I -E, - ikK 



(69) 



where is found self-consistently to satisfy Eq. (1681 . 
After straightforward algebra, one obtains charge neu- 
trality potential in terms of E s , 



V c (k,t) 



Pik • E s (P„ - P)ik • E s 



A- 2 



K(x,i) = - 



P 

4tt 



dV 7 



fc 2 + a;/ 



(70) 



1 % ^V-E s . (71) 



4tt 



The effect of this potential is twofold. Firstly, V c gives 
additional electrical current —<r^W c . In fc-space, 



Ajt(k,t) = -o^ikV c , 
Aj+(k,i) = -a^kV c . 



(72) 
(73) 



Secondly, V c affects diffusion current through Vn^. Af- 
ter some algebra, one obtains 



t 

D 

Alt, 



Aj 



ik 



X 



t 2 







sf 

Ag 



rr 1 



^ o 

4 



k 2 V c 
k 2 V c 



ikVr 



(74) 



which exactly cancels Eqs. (IT2"j) and (j?3")) . Hence, there 
are no additional charge and spin currents attributed to 
V s . 

It is very interesting that one does not need to assume 
any specific form of V c to reach Eq. (|T4"|) . The mathemat- 
ical origin of the exact cancelation is that the additional 
force is conservative (described by ikx (scalar)) and spin 
independent. A spin independent potential generates ad- 
ditional current, but spin accumulation adjusts fast to 
make diffusion current cancel it. Of course, the modi- 
fied spin accumulation n'* affects LLG equation by Eq. 
(|17]). but we claimed that this is ignorable. Consequently, 
any spin independent potential cannot modify the main 
features of our result. As another side remark, this can- 
celation can be verified for exact Coulomb potential even 
without using the charge neutrality approximation Eqs. 

©-CD- 



C. Arbitrary boundary 

Since SMF is usually strongly localized, we have con- 
sidered an infinite boundary problem. However, for a 
finite system, the expression of the charge and spin cur- 
rents should be slightly modified. Note that the key part 
of our theory is to find real space expression of fl^ 1 . (See 
Eq. (|28| ). In real space, fHs a differential operator, 



^real — 



n Q+ n Q 



f i I 



T i I 



-a^v 2 



T i 4- 



Hence, the problem is to find the inverse operator, i.e., 
the Green's function corresponding f2 rea i at the given 
geometry. 

In order to find the Green's function of fl rea i, one 
can get a hint from Eq. (l29l) . Firstly, let Gi(x,x') 
and Gij(x, x') be respectively the Green's function cor- 
responding Laplacian V 2 and modified Helmholtz op- 
erator V 2 — Xjj 2 for the given geometry. In order to 
obtain the Green's function, it is plausible to replace 
1/fc 2 -> -G L (x,x') and l/(fc 2 + Aj/) -> -G ff (x,x') in 
Eq. ([2"9"]>. Then, one obtains 



Gn(x,x') 



A 



^Gl-^Gh ff ; G L +g T G 



(76) 



A 



A 



4 2 



One can verify this is indeed the Green's function of fi rea / 
by showing VtreaiGn = <5(x — x'). Then, one straight- 
forwardly concludes that the generalized expressions of 
charge and spin current [Eqs. (j42 |) - (|43| ] are given by 

j e = P<j-E s -PctV J A'G i (x,x')V'-E s , (77) 

j s = Pie + (1 - P 2 )aE s 

-(l-P 2 )aV J A'G fl (x,x')V'.E s . (78) 

For infinite boundary with vanishing boundary condition, 



47tGl = -l/|x-x'| and 4irG H = -e 
so Eqs. (j42" |) - (j43"l) are reproduced. 



-|> 



*'I/A 3 //|, 



X' 



VI. CONCLUSION 

By constructing spin drift-diffusion equation from the 
equation of motion of conduction electrons, we studied 
the effect of SMF in the presence of spin accumulation, 
spin diffusion and spin flip scattering, which were ignored 
0] or considered only in ID fl9l | in previous theories. 
It turns out that, in realistic regime, the conservative 
part of charge current is canceled by diffusion current, 
and spin current becomes nonlocal. Consequently, the 
magnetization dynamics is affected by nonlocal Gilbert 
damping tensor, instead of the previously reported (lo- 
cal) Gilbert damping tensor. By calculating spin-transfer 
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torque generated by nonlocal spin current, we obtained 
the explicit expressions of the nonlocal Gilbert damping 
tensor. 

Different from the previous work focusing on ID [l9j . 
we also obtained the results for 2D and 3D as well as ID. 
In a ID system, the results of the previous theory are 
reproduced. It turns out that Einstein's relation prevents 
the existence of charge current, but for parameter sets 
which do not satisfy the Einstein's relation, the charge 
current can be induced by the spin motive force. After 
generalizing the results to higher dimension, we find that 
the nonconservative part of SMF plays an important role 
in charge and spin currents. 

As an illustration of suppression of the effect of SMF, 
we demonstrated equations of motion of collective coor- 
dinate of ID current-induced DW motion. We found that 
spin diffusion renormalizes SMF depending on the rela- 
tive magnitude of DW width and spin diffusion length. 

We also investigated the system under spatially ex- 
tended (non-localized) electric field. In this case, it turns 
out that the spatially extended part of the electric field 
can contribute to the charge current even though it is 
conservative. However, the existence of spatially ex- 
tended part of the electric field cannot alter the result 
that the charge current generated by (localized) SMF can- 
not include conservative part. 

Our result is solid in the sense that our principal re- 
sults are not changed by any spin independent potential. 
A spin independent potential modifies current via addi- 
tional electric field, but spin density rapidly adjusts to 
make diffusion current cancel it. Consequently, a spin 
independent potential can modify spin density, but not 
current. 
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where / is the identity matrix. Since 



JW(k,0) = (2n) d / 2 5(k)n 







(A2) 



the second term of Eq. (|A1|) vanishes after applied to 
Af(k, 0). Finally, one obtains e~ nt ^fAf(k, 0) = Af(k, 0). 

Note that the result is obtained without the approxi- 
mation that r s f is small. 



Appendix B: Derivation of Eq. (|27p 



The idea is based on the delta sequence 

lim ne- nt Q(t) = S(t), (Bl) 

n— >oo 

for natural number n. This can be generalized as 

lim ne- uzt S{t) = 5(t), (B2) 

u — > oo 

where z is a complex number satisfying 5R[z] > 0. This 
generalization is obvious in that uze~ uzt is localized at 
t = and uze~ uzt Q(t)dt = 1. 

We now generalize this relation to matrices. For a 
diagonalizable matrix M with eigenvalues Xi satisfying 
Si[Aj] > 0, we claim that 



lim uMe~ uMt e{t) = 5(t)I. 



(B3) 



The proof is simple. Since M is diagonalizable, one can 
write 



uM = Q 



uAi 
u\ 2 



(B4) 
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Appendix A: Absence of Time Variation of 
Equilibrium Number Density 

In this section, we show that e~ nt / Tsf J\f(k, 0) does not 
have time dependence, where fl is given by Eq. (|25l) and 
Af{k,0) = (27r) d / 2 (5(k)(nJ n^) T . It suffices to calculate 
e -nt/r B f a t = because of S(k) factor. It is easy to 
show that fi is idempotent for k = 0. In other words, 
n n \ k = = fi|*=o. Then, for k = 0, 



-nt/r 



Q 7 



n 



t 

n—Q x J 



n=l 



t 



!)• 



(Al) 



for some Q. Then 
uMe~ uMx O(x) 



Q 



\ 

( 5{x) 
8{x) 



i\ 2 e- uX *Q(a 



= S(x)I, 



Q- 1 



(B5) 



as u — > oo. 

By taking M = Q and u = 1/t s /, 

Q 



lim — e - nt/T "e(i) =S(t)I, 



so, for small r. 



T sf 



(B6) 



(B7) 
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which is exactly Eq. (|27l) . 

As a passing remark, one should check that the real 
parts of eigenvalues of are really positive. It is guar- 
anteed by Tr fi > and det O > for any nonzero k. 
However, the approximation Eq. (|27p makes singularity 
at k = 0. Fortunately, this singularity is removed by zk 
factor in Eq. (j2"o) . 



Appendix C: Nonlocal damping tensor and energy 
dissipation 



1. ID damping tensor 



Starting from 



3s 



-\x-x'\/\ sf 

(1-P 2 ) dx' aE s (x',t), (CI) 



2A 



sf 



T dtM-(d x MxM), 



(C2) 



2eM s 3 

spin-transfer torque driven by SMF is given by 



-\x-x'\/\ sf 

x / dx' - d t M' ■ (d x M' x M') 



2A 



M 



x / dx' 



sf 

(1-P 2 ) V ' 



— \x — x | / X s f 



Mf 2A 



sf 



(M x d x M)(M' x ftj,M')] • dtM'. (C3) 
Then, comparing with Eqs. and (f5TT) . 

(1 - P 2 )rj e~\ x ~ x '\/ x 't 
V i:j {x, x') = aSijSix - x') + — 

x(M x d x M)i(M' x d x M')j. (C4) 



given by 

dE 
~dt 



II 



M 



x V ■ dtM 



Y I dx(M x H) • V ■ d t M 

~- J dx(M x H) ■ V • (M x H) 

J J dxdx'(M x H) l 2? y (x,x')(M' x E(%) 



_7 

M, 



up to first order. Here, H' = H(x',t). The first term 
in Eq. (|C4j) gives definitely negative dE/dt. The second 
term, which comes from SMF, gives 

mbm* = Jl^pm ff dxdAu . dxM) 

dt M s JJ V ; 



x^-, (H'.0 x M'). 



2A 



sf 



(C7) 



In order to show dE/dt < 0, it is sufficient to show that 
JJ dxdx' f(x)e~ a][X ~ x 'f(x') is positive for real function 
/. This is verified by Parseval's relation and convolution 
theorem of Fourier transform. 



dxdx'f(x)e- a ^- x '\f(x') 



= / dxf(x)* 



dx'e- a \ x - x '\f{x') 



= J dkF[.f{x)](k)*F J dx'e~ alx - x ' l f(x') 



(A-) 



= V2tt / dkT[f{x)}(k)*T[e- a ^}(k)T[f{x)](k) 



=2aid k mm->o. 



k 2 + a 2 
This implies dE/dt < 0. 



(C8) 



It is interesting to note that 



lim T>ij{x — x) — 5(x — x) 



A sf ->0 



aSi 



(i - p 2 )v 



;(Mx d x M)i{Mx d x M)j\Ch) 



This is exactly the previous result Eq. ([5]) except (1 — P 2 ) 
factor, which should exist regardless of diffusion strength. 

Now, the remaining step is to calculate energy dissipa- 
tion. Energy dissipation is calculated by the integration 
energy density dissipation — H e ff ■ dtM. From now on, 
the subscription e ff is dropped until this section ends. 
Energy dissipation by a nonlocal damping tensor T> is 



2. 3D damping tensor 

From now on, Einstein's convention is used. Compo- 
nentwise expressions of spin current and SMF for a 3D 
system are 



3s,i ®E s i 



P 2 G 

An 



0. I dV^H 



(1-P>, 



4tt 



ft / dV 



|X-X'| 

>-|x-x'|/A a / 



x — X 



E„ 



2eM 3 



d t M ■ (dtM x M). 



—^^•(09) 
(CIO) 
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By integrating by parts in order to remove derivatives in 
front of E, ,-, 



P 2 a 



47T 



(i-p> ; 



4~ 



d t / d A x'E s ^)3' f 



,-|x-x'|/Aa/ 



oE, 



P 2 a 

47T 



47T 



d 3 x'E s , 3 (x')d' i d' J 



-|x-x'|/A„/ 



Here, it is convenient to introduce nonlocal conductivity 
tensor S. 



J dVS y (x,x')£ SJ (x'), 



(C12) 



5 l:) (x,x') = er<5y <5(x - x') 



Air J x — x' 



+ l ; }tT ^—, — . (ci3) 



47T 



X - X' 



Now, spin-transfer torque driven by SMF is given by 

Ms a. • v)m 



eM 



eM t 

M 

m7 



9 4 M / d 3 x' • < S ij (x,x')£ sJ (x') 



(M x cW) 



x y dV • 5y(x,x')(M' x 9,-M') • <9 t M' 



(.C14) 



Hence, the nonlocal damping tensor for 3D is given by 
ZWx,x') = aSijS(x - x') H — j5fc;(x,x') 

x(Mx9 fc M) i (M'xaiM') j . (C15) 

Now, we calculate energy dissipation. This is an ana- 
logue of the previous section. After some algebra, 



dEsMF 

df 



7^7 



d 3 xdVSfcz(x,x') 



x(H-9 fc M)(H'-^M'). (C16) 
In order to show dEsMF / dt < 0, one should verify 

d 3 xd 3 x'{H ■ a 4 M)5y(x,x')(H' ■ djM') > 0. (C17) 



Note that 5y(x, x') is a function of (x — x'). It is con- 
venient to write <Sy = 5y(x — x') at this stage. Similar 
to the previous subsection, by using Parseval's theorem 
and convolution theorem of Fourier transform, 



d 3 xd 3 x'{U ■ d t M)S tj (x - x')(H' • djM') 



= / d 3 fc^[H-a 4 M](k)*5y(k)J-[H-a 3 M](k()C18) 



Here, Sij (k) is given by 



(2tt) 3 / 2 



k 2 



k2 + Kn 



(C19) 

Note that the integrand in Eq. (IC18|) is an expecta- 
tion value of 3 x 3 matrix Sij (k) with respect to vector 
J-[H ■ 9jM](k). Since Sij(k) is a Hermitian matrix, the 
integrand is always positive if the eigenvalues of <Sy(k) 
are positive. Recall that the eigenvalues of matrix kikj 
is given by (0, 0, k 2 ). The corresponding eigenvectors are 
definitely the eigenvectors of Sij(k), so the eigenvalues 



of 5y(k) is given by yjj^ , j^^, (2ff)3/a fc2+A _. 
which are all positive. This proves dE/dt < 0. 



Appendix D: ID DW motion in the presence of spin 
diffusion 



1. Collective coordinate equation of ID DW for 
space-time dependent spin current and spin density 



The starting equation is ID version of Eq. (|47|) . 

1 71cUr\ <9M 
1 + ' 1 

1 



1 + /3 2 M s y ^ = -7oM x H e// 



a l + P 2 M s J M~ M X dt 
+b dM ^ M ^ <9M 

9a; M a 5a; ' 



(Dl) 



where bj = ^iBjs/sM s . The main difference from the 
conventional approach is that bj and n s are space-time 
dependent. The equation is rewritten as 



7o M x Htot = 0, 



(D2) 



where 



7 H tot = - 1 



1 



n s (J,B 



1 + /3 2 M s J M| ' ' dt 



M dM 

— x —+ 7oH e// 



1 + f3 2 M s J M s dt 
M dM /3bj dM 



+hj W 2 X dx 



Ms dx 



(D3) 



We set the effective magnetic field as H e ff — jjsijjsr + 

^-jf^x + H^, where A and represents exchange cou- 
pling and anisotropy, and is dipole field. In addition, 
the magnetization profile is set to be tail-to-tail trans- 
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verse wall, 



M = M s (cos 9, sin 9 cos tj>, sin 9 sine 

'x-x{ty 



sin 9 = sech 



cos# 



tanh 



X 

x-X(t) 
A 



(D4) 
(D5) 

(D6) 
(D7) 



where X(t) is DW position and 0(t) is tilting angle. 

Note that Eq. (|D2p implies H tot = aM for some a. 
Then, one can define force density / = — H tot • and 
torque density t = — Htot ■ 75^, which are identically 
zero. Then, collective coordinate equation is given by 
calculating total force and total torque, 



F 



fdx = 0, 



T = rdx = 0. 



(D8) 
(D9) 



Each equation implies respectively, 

d<p adX 1 / d<f> (3dX 
~~dt + X~dt + 1 + /3 2 \~dt + X~dT 

J M s \ X J 2X 

= -^--xJ bjSech [— J2V (D10) 
i ax 9^ i /iay 90 
a at + a at + 1 + p 2 1 a at + " at 



sech 2 

AT, 



x — X \ 
A J 2A 



A 



A 



A 



2A M, 



sin(2$.l) 



where 6j is space-time independent part of bj which 
comes from applied spin current, bj = bj — bj, Kd cor- 
responds integration of dipole field. Here, a and 7 are 
renormalized parameter by the same way with nj^. One 
can check that Eqs. (|D10[) and (|D11[) reproduces Eqs. 
([531 and ([51]) if An s = and 6} = 0. 



2. ID DW motion in the presence of SMF and spin 
diffusion 



Now, we apply 

bS = ^ B a 1 - P 2 
3 eM s 2X sf 



dx'e-\ x - x 'V x «f E s (x' ,f](pi2) 



h , 2 (x-X\d<f> 



An, 



ar s/ 1 - P 2 k 2 + P„PA s/ 2 t fcff s 
eA^l-P 2 P + A;/ P~' 



(D13) 
(D14) 



to Eqs. (|T3T0|) and (|DTT]) . 

First, we calculate the integral J An ,5 sech 2 ( a T y ) dx 
which corresponds to the effect of spin density . Note that 
Au s is an odd function of x — X. Hence, the integrand 
is an odd function so the integral vanishes. 

The next step is to calculate / sech 
ter some algebra, 



2 ( x-X \ dx 
\ A I 2A 



M. Af- 



6 s 

-^sech' 
A 



x — X\ dx 



2A 



A 



(l-P>C g0 [[ e -CW-u'\ sech 2 usech 2 u , dudu , 

4A 2 at 77 
(i-p 2 )7?ca0 f 7rCfc 2 csch 2 



4A 2 at 
(1 - p 2 )t/C a0 ' v 



fc 2 + c 2 



27rCfc 2 csch 2 ^' 

4 a- at./ PTc 



£7T \ 

2 -±dk 



(D15) 



where £ = A/A s y, u = (2; — X)/X, and u' = (x' — X)/X. 
Parseval's relation and convolution theorem of Fourier 
transform is used at the third step. By using the identity 



^ — nnk 



(Die) 



the integral can be expressed by Laplace transform C 



2ttCP csch 2 
k 2 + ( 2 



dk = J2 C 



= E 



iir^nk 2 



(mr) 



k 2 + c 2 

16 sinmr^ 



^T[Jo wr(u+l) 
Recall the Fourier series of a sawtooth function 
16 



m.7) 



— sin nwx = 8(2n + 1 — x) 



(D18) 



for 2n < x < 2n + 2. Then the integral is given by 

16 sinwrCit 

-du 



y , 

, Jo n7r ( u + 1)' 

n— 1 u v ' 



E 



— 8(2n+ 1 - uQ 



n=0 " T 



du 



1 



(D19) 



By using the relation 
d n T'(z) 



dz n r(z) 



(-1)" +1 «!E 



(D20) 
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one obtains the integral in closed form as 
8 F(C) _ [°° 27rCfc 2 csch 2 (*£ 



3 C 



= -4 - 4C + 2( 2 



k 2 +( 2 
T(z) 



dk 



-4 - 4C + 2C 



r"[§)r 



(0 



.(D21) 



The prefactor 8/3£ is to make liniA 3/ ^o ^(C) = 1- There- 
fore, we finally obtain 



sech : 

A 



x — X\ dx 
2A J 2A 



(D22) 



and, consequently, 



i ax 



a at 



at 



~A~ ~ 3A2 ~di: 

2f) d<j> jK d 

A ~" " 



(D23) 



where 77 is renormalized parameter by spin diffusion de- 
fined as fj = (1 - P 2 )F(()ri. 
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